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Dynamics of spin-polarized neutral particles, such as neutrons or neutral atoms and molecules, 
in magnetic fields is studied. A new regime of motion is found where particles move mainly in one 
' direction forming a well-collimated beam. This regime suggests a mechanism for creating devices 

emitting directed beams of neutral particles. 
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The procedure of creating and accelerating directed beams of particles is of high practical importance. It would be 
difficult to enumerate various applications of such directed beams. Charged particles, as is known, are well manipulated 
by electric and magnetic fields [1]. To the contrary, for neutral particles there are no known ways of creating and 
accelerating directed beams by means of nonresonant electromagnetic fields. Directed neutral beams have been formed 
by employing mechanical collimators selecting particles from an isotropic distribution, by accelerating neutral particles 
through long tubes with a high pressure difference between the ends, as is done in molecular beam masers [2-4] , and 
for resonance atoms by using laser beams. In this paper we show that there exists a magnetic mechanism that allows 
the formation of directed beams of neutral particles with strong acceleration. This mechanism requires no mechanical 
collimators, but does require a particular configuration of magnetic fields and a special polarization of particle spins 
at the initial time. A brief announcement of the results presented below has been given in Ref. [5]. 

To demonstrate the principal possibility of this new dynamical mechanism, we consider a rarefied gas of particles 
whose collisions at the first approximation can be neglected. Let particles be placed in magnetic fields whose space 
variation is sufficiently smooth, so that the semiclassical description can be applied [6]. Then for a particle with 
mass m and magnetic moment /i , one may write the evolution equations for quantum-mechanical averages. For the 

average real-space variable R= {R x , R y , R z } one has 

d 2 R a Mo - SB 

with initial conditions R (0) —Ro and R(0) = Rq, where the dot means, as usual, a time derivative. It is easy to check 
[6] that Eq. (1) follows from the definition of R as an average of a position operator with a wave function satisfying 
the time-dependent Schrodinger equation. The average spin S= {S x , S v , S z } satisfies the equation 
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d S Mo -* -> /„n 

with an initial condition 5 (0) = {Sq , <Sq , So}- Let us take the total magnetic field B=Bi + B2 as a sum of two 
terms, 



B 



= Bj (V l x +R y l y +\R Z e z ) , B 2 = B 2 h (i), (3) 



the first being a quadrupole field, parametrized by its gradient B[ and the anisotropy parameter A, and the second, 
a transverse field, with 

h (t) = h x ~e x +h y ~e y , I /i |=1, (4) 

depending only on time but not on real-space variables. Such magnetic fields are easy to form and arc often used 
in different applications. For example, the quadrupole fields are the basis of quadrupole magnetic traps and the 
Ioffe-Pritchard traps [7,8]. A transverse rotating field, with h x = cosut and h y = smut, has also been employed in 
magnetic traps [9]. This rotating field will be used below as a concrete example. However, the effect we consider 
exists not solely for such a rotating field but for a wide class of transverse fields satisfying the condition 

d ~h /dt < Li B 2 /h . (5) 

To be absolutely concrete, we also take the anisotropy parameter A = — 2 so that V • B\= 0. For what follows, it is 

convenient to measure the components of the space vector R in units of the characteristic length R = B 2 /B[ defining 

the radius of the field zero in the radial direction. To this end, we define the dimensionless vector r=R / R = {x,y,z}. 
Introduce the characteristic frequencies 

lu\ = /j, B' 1 /mR , lo 2 = /j, B 2 /h, (6) 

the first of which, w\, corresponds to the motion of particles in real space and the second, oj 2 , to the spin motion. 
The physical meaning of these frequencies becomes evident after we present Eqs. (1) and (2) in the form 

= J* (> e x +Sy t y -2S Z e 2 ) , ^ = u 2 A S, (7) 

where the matrix A = [A a p], with a,(3 — 1,2,3, consists of the elements A aa = 0, A i2 = —A 2X = — 2z, A i3 = 
—A31 = —y—hy, A 2 3 = —As2 = x + h x . From (7) it is really clear that co\ and ui 2 are the characteristic frequencies 
of the space and spin motions, respectively. 

Notice that the system of equations (7), with a nonuniform magnetic field, is invariant under the change S— * — S 
and r — > — r . This invariancc can be called the Stern-Gerlach symmetry since in the particular case of a uniform 
magnetic field one would recover the conditions of the Stern-Gerlach experiment. 

To solve the system of nonlinear differential equations (7), let us recall that (1) and (2) are derived in the semiclassical 
approximation whose criterion of validity is the slow space variation of magnetic fields [6] . In our notation, this criterion 

can be written as the inequality |u>i/w2| <§C 1. This condition shows that the space variable r can be treated as slow, 

compared to the fast spin variable S- From inequality (5) it follows that h is also slow as compared to S- For the 
rotating field, condition (5) simply means that w«W2. Therefore, the system of nonlinear equations (7) can be solved 
by employing the method of scale separation [10-13] which is a variant of the Krylov-Bogolubov averaging method 
[14,15]. A detailed description of this approach as applied to the nonadiabatic dynamics of atoms in nonuniform 
magnetic fields has been given in Refs. [16,17]. Following the method of scale separation, we, first, need to solve the 
equation for the fast variable, keeping there the slow variables as quasi-integrals of the motion. Then, under fixed r 
and h, the second equation from (7) can be solved exactly. The resulting solution is 
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S(t) = J2 Oi Si (t), a t =S (0)- (0), S t (t) =b t (t) exp {&(*)} , (8) 



b l (t) 



(a 4 Ai 3 + A 12 A 23 ) ^ x + ( ai A 23 - A 12 A 13 ) 7 y + (of + A\ 2 ) e . 
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Ci = (|«,| 2 - A\ 2 Y + (M 2 + A\ 2 ) (Al, + Aj 3 ) , c? = A\ 2 + A% + Aj 3 , 

ai )2 =±ia, a 3 = 0, /%(*)= w 2 / a^t) dt , /? 3 (i) = 0. 

Jo 

The fast solution (8) is to be substituted into the equation for the slow variable, averaging the right-hand side of it 

over an interval of time much longer than the period of fast oscillations 27t/cj 2 , which gives (S) = 03(63), and for the 
case of the rotating field 

_ (1 + x)S§ + yS^ - 2zS§ -* _ (x + cosa;t) ~e x +(y + smut) ~e y ~2z e 2 
° 3 ~ [(1 + x) 2 + y 2 + Az 2 } 1 / 2 ' b3 ~ [1 + 2(x cos cot + ysincot) + x 2 + y 2 + Az 2 } 1 / 2 ' U 

In this way, we obtain 

, _ ^ _ x w? K 1 + a; ) 5 'o + vSo - 2z Sq] (x e x +y e y +8z e z ) 

{ul[S x e x +Sy e v -2S z e z )) = ^ — (10) 

2 [(1 + 2x + x 2 + y 2 + Az 2 ){\ + x 2 + y 2 + Az 2 )] 1 ' 2 



As a result, we come to the equation 

d 2 7 



dt 2 



--< wl(S x e x +S y 7 y -2S Z e z ) > (11) 



describing the averaged motion of particles. If we take adiabatic initial conditions for the spin polarization such that 
Sq 7^ and Sq — S§ — 0, then for r<c 1, the right-hand side of Eq. (10) reduces to the simple harmonic force 
\lo 2 Sq(x e x +y e y +8z e z ). Hence Eq. (11) describes the adiabatic motion of particles in a harmonic potential. For 
< 0, this is the standard harmonic oscillations of trapped particles, while for Sq > 0, the latter are not confined 
and escape from the trap in all directions. Neither of these known cases is of interest for us. Our aim is to find a 
principally different regime of motion, when the particles are neither completely trapped nor uniformly coupled out 
but move preferably in one direction forming a collimated beam. Assume that at t — the particles are prepared in 
the state with an initial spin polarization 

S x q - St = 0, S z q = -S , (12) 

which is referred to nonadiabatic initial conditions [18]. The ways of preparing such polarized initial states are 
discussed in detail in Ref. [18]. To take into account the finite size of a device, we introduce the device form-factor, 
tp(r), and use the notation 

f( 7) = ri 3 ) 

Jy ' [(l + 2x + x 2 +y 2 + Az 2 )(l+x 2 +y 2 +Az 2 )} 1 / 2 ' K ' 

Then the motion of particles is described by the equation 

Suj\fz (x e x +y e y +8z e z j . (14) 



dt 2 



From here it is seen that the motion along the x and y axes is similar to each other. One may also notice that Eq. 
(14) is invariant under the inversion z — > — z and S — > — S. Thence, the trajectories for S > are mirror-symmetric, 
with respect to the x — y plane, to the trajectories for S < 0. Therefore, it is sufficient to study only a case, say, when 
S > 0, which will be assumed in what follows. It is also tempting to simplify the equations by passing to dimcnsionlcss 
time measured in units of (y/SuJi)^ 1 . For this purpose, we define r = VSuit. Thus, from (14), we come to the system 
of equations 

% = s ^ a?- 8 '* (15) 

An analytic solution of Eq. (15) is possible only for the initial stage of the motion, when | r | <C 1. The 
corresponding solutions are given by the Weierstrass and Lame functions [19,20]. From the properties of these 
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functions, it follows that the axial motion of particles is bounded from below by the minimal value z m i n = (zq — C) 1 ^ 3 , 
with C = (3/16)zq, and zo = z(0) being an initial axial position. So, particles can escape only in the positive z— 
direction. The motion along the axial direction is much faster than along the radial direction. This means that 
beam collimation begins already at the initial stage of the process. To analyze accurately the whole motion, we 
have to resort to numerical calculations. Let us specify the form-factor tp( r ) assuming a spherical device with 
ip(r ) = cxp(— | r \ 2 /L 2 ) , in which L is a characteristic size of the device, e.g., the radius of a coil forming magnetic 

1 /2 

fields. The maximal velocity in the axial direction is w max = (io + . Therefore, taking L sufficiently large, 

it is feasible to get arbitrary strong acceleration in the axial direction. At the same time, the maximal velocity in the 
radial direction, in the case of L 3> 1, is much less than w max , which can be accepted as the definition of collimation. 

The results of numerical calculations are presented in characteristic figures corresponding to the initial position at 
the center of the device and to initial velocities varying in the interval [—1,1] in all directions. All figures for different 
L>1 are qualitatively the same, because of which we fix here L = 1000. Fig. 1 shows the trajectories of particles in 
the x — z plane at the beginning of motion. It is clearly seen how the particles having initially negative velocities in 
the axial direction, after reaching the minimal value z m i n , turn to the positive z direction. In Fig. 2 the trajectories 
are shown for longer times, demonstrating how a well-collimated narrow beam is formed, being stretched in the axial 
direction more than an order of magnitude stronger than in the radial one. In Fig. 3 the velocities v(t) = x(t) and 
w(t) = z(t) are pictured, illustrating the acceleration process in the axial direction. When varying L, all figures 
remain qualitatively the same. The sole thing that changes is the scale on the z-axis. Increasing L by an order 
squeezes the z scale approximately twice in the figures with trajectories and three times in the figures with velocities. 
Thus, the maximal velocity in the axial direction in the dimensionless units is w max = 6 for L = 10; w max — 19 for 
L = 100; and w ma x — 60 for L = 1000, in accordance with the law w ma x — {2^/ttL) 1 / 2 . 

We would like to recall that equations (15) are derived from (7) which is identical to the initial Eqs. (1) and (2). 
In the derivation of (15) from (7) the sole assumption used is the existence of the small parameters |wi/w 2 | <C 1 
and \oj/uj2\ <SC 1. This made it possible to apply the scale separation approach [10-13,16,17] whose mathematical 
foundation is based on the Krylov-Bogolubov averaging method [14,15]. 

In order to understand how to choose the characteristic parameters B[ and B 2 of the magnetic field (3), it is 
necessary to return to dimensional units. For simplicity, we set S ~ 1. Then, the maximal z velocity is w max — 
(2y/nLfioB2/m) 1 / 2 . To achieve an effective acceleration for a given sort of particles with fixed no and m, we should 
take L > 1 and sufficiently large B 2 . However, increasing L and B 2 implies the increase of the size of a device 
I = LB 2 /B' l . Hence, to achieve L » 1 for a given device requires that B 2 <C IB[. At the same time, the existence 
of the small parameter (wi/o^l -C 1 yields h~ 2 (B[) 2 / 'm[i B 2 <C 1. Therefore, for a given sort of particles and a given 
device we must have {% /m/xo) 1 ^ 3 (-Bi) 2 ^ 3 < B2 ^ IB[. This fixes the required relation between the parameters of 
the magnetic fields. The frequency of the rotating field, according to Eq. (5), is to be much smaller than tu 2 . 

The effect described in this paper is rather general, and choosing the corresponding parameters, one could realize 
such a directed acceleration for any kind of neutral particles having spins. To show that the parameters to be chosen 
are quite realistic, let us make numerical estimates for 87 Rb in a magnetic trap of Ref. [9]. Then, the gradient of the 
quadrupole field is B[ = 120 G/cm and the amplitude of the rotating field is B 2 = 10 G, with the rotating frequency 
co = 5 x 10 4 s _1 . This gives uj\ <~ 10 2 s~ 1 and ii)2~5x 10 7 s _1 . Hence, the inequalities wi « w « W2 hold true. The 
required small parameters uj\juj 2 ~ 10~ 6 and lo/uj 2 ~ 10~ 3 are really very small, because of which the scale separation 
approach provides very accurate solutions differing from the exact ones by negligible corrections. The radius of the 
atom cloud is R <~ 0.1 cm. Taking L = 10, the characteristic radius of coils would be / <~ 1 cm. For L = 100, this 
would be I ~ 10 cm. The maximal velocity w max w 60 cm/s for L = 10 and w max w 200 cm/s for L = 100. Starting 
from an isotropic distribution of velocities, io ~ io, as a result of the preferable acceleration along the z axis, one can 
obtain the radially squeezed velocity diagram with a rather large squeezing factor. As the phase portraits in Figs. 2 
and 3 show, the squeezing factor is about 20 for velocities and 50 for the real phase variables. Let us stress that such 
a high degree of squeezing is achieved under a rather unfavorable assumption of a spherical device with a spherically 
symmetric form-factor <p(r). Taking a cigar-shaped device would strongly enhance the degree of collimation. It is 
also possible to achieve additional squeezing of the beam taking a quadrupole field in Eq. (3) with essentially different 
field gradients along the axial and radial directions. Thus, one can reach quite high degree of collimation with a 
squeezing factor of 100 or 1000 and more. The realization of the effect we described is quite feasible by using existing 
magnetic traps. What one needs to do is to prepare particles in the initial state with polarization (12). This could 
be achieved in several ways. For example, one could prepare particles in the desired state inside one trap and then 
quickly load them into another trap with the considered magnetic field [18]. Or it might be possible to form the 
necessary initial polarization by means of a short pulse. 

The presented calculations are based on single-particle trajectories. Many-particle physics, related to collisions of 
particles, has not been treated. Since we would like to stress the practical applicability of the considered mechanism, 
it is useful to give at least a coarse estimate on what happens if collisions are included. And the most important is to 
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define conditions when the collisions of particles do not essentially disturb the single-particle picture and permit one 
to realize the semiconfining regime of motion, when particles move predominantly in one direction. 

In order to include particle collisions, we need to add to the right-hand side of the evolution equation (14) an 
additional term describing an effective force caused by these collisions. It is customary to treat the collisional force 

as a random variable. To this end, we can model this force by a random vector 7 £ , where 7 is a collision rate and 

£ = {^x,^y,^z} is a stochastic vector variable. Following again the common way, we may interpret the set {^(i)}, 
with fi = x, y, z, as a set of Gaussian random variables characterized by the stochastic averages 

« »= , « (i 7 ) »= 2D^6^6(t - t') , (16) 

in which is a diffusion rate in the fx direction. Adding the random collisional force to the right-hand side of the 
evolution equation (14), we have, instead of Eq. (15), the system of equations 

d 2 x 

— = Sujjfzx + 7^ , 
d 2 y 



^-=Suj 2 fzy + rt y , (17) 
= 8Su; 2 fz 2 + rt z . 



dt 2 

As is evident, if particle collisions are intensive, so that the motion is dominated by the random collision terms, 
then no organized motion of particles coherently moving in one direction is possible. The semiconfining regime can 
survive only if influence of collisions is sufficiently weak, so that the random terms in Eqs. (17) can be treated as 
perturbation. In such a case, the solutions to Eqs. (17) can be presented as 

x = xi + x 2 , y = yi + V2, z = zi + z 2 , 

where xi, yi, and z\ are the solutions of the unperturbed Eqs. (15) and x 2 , y 2 , and z 2 are the solutions to the 
equations 

' 2 = Sujjf (Z1X2 + X1Z2) + lix , 



dt 2 

d 2 y 2 
dt 2 



= ^ 2 /(ziy 2 + 2 /iz 2 )+ 7 ^ , (18) 
= l6Su;lfziz 2 +7^ . 

The unperturbed functions x\, y\, and z\ can be considered as slow compared to the random functions x 2 , y 2l and 
z 2 . Keeping in Eqs. (18) the unperturbed functions as quasi-invariants, one obtains the solutions 



x 2 (t) = f G x {t - t') [jUt') + Soj 2 fx lZ2 {t')\ dt' 
Jo 



z 2 (t) 



I G z {t-t')lUt')dt' 
Jo 



where the solution y 2 (t), having the form similar to x 2 (t), is not written down and 

^ , , smh(u) x t) „ , . sinh(w z t) fZ 9 ' 

G x (t) = G z (t) = LJLi, u x = JSuj 2 fz 1 , lo z =4lj x 

LU X UJ Z v 

According to the stochastic averages (16), we have 

< x 2 (t) >=< z 2 {t) >= , 
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< x\{t) »= 



l 2 D x t 



sinh(2w a ,t) 
2u> x t 



+ 



wt.x 2 ^ 2 D z t f smh.(w z t) 
H — oxo 9 i cosh(w x t)cosh(w z i) H — [cosh(w z t) - cosh(u; x i)J 



u z t 



— — smh(uj x t)sm]i(uj z t) — If, 



(19) 



< z; >= 



7^ 



sinh(2u; z i) 



- 1 



2to z t 

As is seen from here, the collisions will not disturb much the ordered motion of particles provided that 

7 2 £ 



< 1 , D = sap{D x , D y , D z } 



If we take for estimates the collision rate as 7 ~ hpa /m, where p is the density of particles and a is a scattering 
length, and the diffusion rate as D <~ hsT/fr, where T is temperature, then we get the inequality 

^f- « 1 . (20) 

The latter shows that the influence of random particle collisions, disturbing the organized scmiconfincd motion, can be 
negligible if density, temperature, or the scattering length are small enough to satisfy condition (20). When inequality 
(20) does not hold, the organized directed motion of particles will be essentially spoiled by collisions. Then the motion 
becomes more complicated, at the same time becoming of no interest for our purpose. Our aim here has been to find 
the conditions when the directed semiconfined motion of particles is possible. Such a regime looks like feasible since 
one always can satisfy condition (20) by varying the parameters of the system. 

Concluding, we have advanced a novel general mechanism for creating well-collimated beams of neutral particles 
by means of magnetic fields. Such particles could be neutrons or neutral atoms and molecules with nonzero spin. In 
particular, such a mechanism can be employed for creating narrow beams of molecules in molecular-beam masers or 
directed beams of neutral atoms for other purposes. The mechanism does not depend on statistics and can be used 
for Bose as well as for Fermi particles. Varying the magnetic field parameters and the shape of a device, one can 
regulate the beam characteristics in wide limits, achieving the desired degree of collimation. 
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Figure Captions 

Fig. 1. The trajectories of particles at the initial stage of the acceleration process, for < t < 5. 
Fig. 2. The trajectories of particles for < t < 20. 

Fig. 3. The velocities of particles in the radial, v(t), and axial, w(t), directions for < t < 100. 
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